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Abstract 

A review of the superstatistics concept is provided, including var- 
ious recent applications to complex systems. 

1 Introduction 

Complex systems often exhibit a dynamics that can be regarded as a 
superposition of several dynamics on different time scales. As a very 
simple example consider e.g. a Brownian particle moving through a 
changing environment. Assume that the environment exhibits tem- 
perature fluctuations on a large scale. Then there is a relatively fast 
dynamics given by the velocity of the Brownian particle and a slow one 
given by the temperature changes of the environment, which is spatio- 
temporally inhomogeneous. The two effects produce a superposition 
of two statistics, or in a short, a 'superstatistics' [Il[2l[3llll[5l[6l[7l[8]. 
The concept of a superstatistics was introduced in P], in the mean time 
many applications for a variety of complex systems have been pointed 
out piillllliaillllllllSlligiiailHlllSlsiE^ The sta- 

tionary probability distributions of superstatistical systems typically 
exhibit non-Gaussian behaviour with fat tails, which can decay e.g. 
with a power law, or as a stretched exponential, or in an even more 
complicated way. 

Essential for the superstatistical approach is the existence of an in- 
tensive parameter (3 which fluctuates on a large spatio-temporal scale. 
For the above simple example of a superstatistical Brownian parti- 
cle, (5 is the fluctuating inverse temperature of the environment, but 
in general (3 can also be an effective friction constant, a changing 
mass parameter, a changing amplitude of Gaussian white noise, the 
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fluctuating energy dissipation in turbulent flows, a fluctuating volatil- 
ity in finance, an environmental parameter for biological systems, or 
simply a local variance parameter extracted from a signal. Some su- 
perstatistical models exhibit anomalous transport, others don't. Most 
superstatistical models are somewhat 'less anomalous' than Levy-type 
models, in the sense that usually more of the higher moments exist 
as compared to Levy processes. The tails of the distributions exhibit 
'fat' tails, but usually these are less pronounced than for a Levy dis- 
tribution. 

The superstatistics concept is very general and has been applied to 
a variety of complex systems. Recent successful applications include 
hydrodynamic turbulence [H [101 El [HI [12] , pattern forming systems 
[Hj, cosmic rays [15], solar flares [13], mathematical finance [T7 1 [T8 l [T9] . 
random matrix theory [20], networks [21], quantum systems at low 
temperatures [23], wind velocity fluctuations [13125], hydro-climatic 
fluctuations p2] and delay statistics in traffic models [24]. The aim in 
the following is to explain the basic concepts and applications in an 
easy-going way. 

2 The basic idea 

Consider a complex system in a stationary nonequilibrium state that 
is driven by some external forces. Usually we think here of a physical 
system (e.g. a turbulent flow) but we may easily apply similar tech- 
niques to economic, biological, social systems, where the meaning of 
the mathematical variables will be different, though the mathematical 
structure is similar. Generally, a complex system will be inhomoge- 
neous in space and in time. Effectively, it may consist of many spatial 
cells (or, the measured time series may consist of many time slices) 
where there are different values of some relevant system parameter f3. 
The cell size is effectively determined by the correlation length of the 
continuously varying /3-field. Superstatistical systems are character- 
ized by a simplifying effect, namely the fact that the relaxation time 
is short so that each cell can be assumed to be in local equilibrium 
(in a certain approximation at least). Sometimes this property will be 
satisfied for a given complex system, sometimes not. 

In the long-term run, the stationary distributions of the supersta- 
tistical inhomogeneous system arise as a superposition of Boltzmann 
factors e~^^ (or analogues of Boltzmann factors describing the local 
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behaviour of the system under consideration) weighted with the prob- 
abihty density /(/?) to observe some value /? in a randomly chosen 
cell: 

r-oo 1 

p{E) = f(p)-^^p(E)e-f'^dP (1) 

Here E is an effective energy for each cell, Z(P) is the normalization 
constant of p{E)e~^^ for a given (5, and p{E) is the density of states. 

A typical example is a Brownian particle of mass m moving through 
a changing environment in d dimensions. Such a Langevin model has 
applications in many different areas of science. In the simplest case 
we may write down a linear local Langevin equation for the velocity v 

ij = -jv + aL{t) (2) 

(L(t): d-dimensional Gaussian white noise) which becomes supersta- 
tistical due to the fact that for a fluctuating environment the parame- 
ter 3 := becomes a random varaible as well: It varies from cell to 
cell on a rather large spatio-temporal scale T (this time scale T should 
not be confused with the temperature, which is denoted as /3~^). Of 
course, for this example E = ^mv"^, and while the local stationary 
distribution in each cell is Gaussian 

Pirn = (J^Y' e-h^"'^\ (3) 

the marginal distribution describing the long-time behaviour of the 
particle 

P{v) = / f{f3)p{mdP (4) 



Jo 

exhibits non-trivial behaviour. The large- tails of this distribution 
depend on the behaviour of /(/?) for /3 ^ [1]. As a result of the 
integration over /3, the probability distribution p{v) will typically have 
fat tails. These can be e.g. power law tails, stretched exponentials, or 
whatever. 

One of the most important example for practical application is a 
X^-distribution for /(/?) (see next section for concrete formulas) and a 
Gaussian distribution for p{v\P). In this case one obtains from eq. ^ 

P{v) ^ (5) 

(1 + (g- l)i6t;2)'?-i 

where q and b are suitable parameters, and v = \v\. The function on 
the right-hand side of eq. ([5]) is called a g-Gaussian [34] and denoted 
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hy Cq ^ . Note that g-Gaussians decay asymptotically with a power 
law for V ^ oo. They reduce to ordinary Gaussians for q ^ 1. But 
power laws are not the only possibility one can get for superstatistical 
systems. Many examples will be discussed in the following sections. 

A generalized thermodynamics for superstatistical systems has been 
recently developed in [8] (see also [5] for early attempts). Here one 
starts from a Boltzmann-Gibbs-Shannon entropy function that also 
includes contributions from the fluctuations in /3. Given that entropy, 
one can do formally thermodynamics. Ordinary thermodynamics (no 
fluctuations in /?) is contained as a special case in this more general 
formalism. 

3 Typical distributions f(f3) 

The distribution /(/?) is determined by the spatio-temporal dynamics 
of the driven nonequilibrium system under consideration. By con- 
struction, P is positive, so /(/?) cannot be a Gaussian. Let us here 
consider important examples of what to expect in typical experimental 
situations for driven nonequilibrium systems. 

a) There may be many (nearly) independent microscopic random 
variables variables S^j, j = 1, . . . , J contributing to /3 in an additive 
way. For large J their rescaled sum Z]/=i ^^^^ approach a Gaus- 
sian random variable Xi due to the Central Limit Theorem (CLT). 
In total, there can be many subsystems consisting of such microscopic 
random variables, leading to n Gaussian random variables Xi, . . . , X„ 
due to various degrees of freedom in the system. As mentioned be- 
fore, (3 needs to be positive and a positive /3 is obtained by squaring 
these Gaussian random variables. The resulting /3 = "^^i Xf is x^- 
distributed with degree n, i.e. 

The marginal distributions obtained by integrating over all /? exhibit 
power-law tails for large enerfies E. They are g-exponentials, p{E) ~ 
e~^^ = {1 + {q — l)bE)^^/^''"^\ where q and b can be related to n 
and Pq [261 I27j . Note that this statistics arises as a universal limit 
dynamics, i.e. the details of the microscopic random variables (e.g. 
their probability densities) are irrelevant. 
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b) The same consideration as above may apply to the 'temperature' 
rather than (3 itself. may the sum of several squared Gaussian 

random variables arising out of many microscopic degrees of freedom 
^j. The resulting is the inverse ^^-distribution given by 

It generates distributions that have exponential decays in [28^ 
m [29]. Again this superstatistics is universal: details of the are 
irrelevant. 

c) Instead of (3 being a sum of many contributions, for other sys- 
tems (in particular turbulent ones) the random variable /5 may be 
generated by multiplicative random processes. We may have a local 
cascade random variable Xi = n/=iO' where J is the number of 
cascade steps and the are positive microscopic random variables. 
Due to the CLT, logXi = Ylj=i^^S^j becomes Gaussian for large J 
if it is properly rescaled. Hence Xi is lognormally distributed. In 
general there may be n such product contributions to the supersta- 
tistical varaiable f3, (3 = HiLi -^i- Then log/3 = Y17=i ^^S^i is a sum 
of Gaussian random variables, hence it is Gaussian as well. Thus /? is 
lognormally distributed, i.e. 

/(/?) = ^exp[-c(ln/?-6)2], (8) 

where a, b, c are suitable constants. The result is independent of the 
details of the microscopic cascade random variables , hence there is 
universality again. This type of lognormal superstatistics is particu- 
larly relevant for turbulent flows [12} [TOl [911301 [TT] . 

4 Asymptotic behaviour for large en- 
ergies 

Superstatistical probability densities, as given by eq.([T]) or typi- 
cally exhibit 'fat tails' for large E, but what is the precise functional 
form of this large energy behaviour? The answer depends on the distri- 
bution /(/?) and can be obtained from a variational principle. Details 
are described in [4], here we just summarize some results. For simplic- 
ity, let us put p{E) = 1 in eq. ([T]). We may define a new probability 
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density / by 

m:-cm, (9) 

where c is a suitable normahzation constant. The new density / ab- 
sorbes the /3-dependence of the local partition function Z(P). With 
this notation, p{E) can now be regarded as the Laplace transform of 
/. Renaming / — > / we obtain 

/•oo 

p{E) ~ / me-f^dp 







oo 



g-/3£+ln/(/3)^^ 



^ gSup^{-/3£;+ln/{/3)} 
^ g-feS+ln/Cfe) 

= /(/3i^)e-^-^. (10) 

Here we used the saddle point approximation. /?£; is the value of /3 
where the function —(3E + ln/(/3) has a maximum. The expression 

sup{-/3i? + In /(/?)} (11) 

corresponds to a Legendre transform of ln/(/3). 

For the case where /(/?) is smooth and has only a single maximum 
we can obtain the supremum by differentiating, i.e. 

snp{-f3E + \nf{(3)} = -PeE + lnf{/3E) (12) 

where Pe satisfies the differential equation 

= -E + {lnfi(3)y = -E+j^. (13) 

By taking into account the next-order contributions around the max- 
imum, eq. (fTO|) can be improved to 

,(E) ~ fML^^L. (14) 

/(&))» 

Let us consider a few examples. Consider an /(/3) which for small 
P is of the power-law form /(/3) ~ P^, 7 > 0. An example is the x^- 
distribution of n degrees of freedom, which was mentioned previously: 



(/9o > 0, n > 1). This behaves for /3 ^ as 

~ (16) 

i.e. 

7 = ^-1. (17) 

Other examples exhibiting this power-law form are the so-called F- 
distributions[Tl[29]. With the above formalism one obtains from eq. (jl3|l 

Pe = 1 (18) 

and 

p{E) ~ E-"'-\ (19) 

These types of /(/?) form the basis for power-law generalized Boltz- 
mann factors (q-exponentials) used in generalized versions of statis- 
tical mechanics, so-called non-extensive statistical mechanics |31l l32| 
[33} . These depend on an entropic index q coming from a more 
general entropy functional, and the relation between 7 and q is 

7 + 1 = ^. (20) 
q-1 

Another example would be an /(/5) which for small /3 behaves as 
/(/3) ~ e~'^/^, c > 0. In this case one obtains 

& = ^ (21) 

and 

p{E) ~ £;-3/4g-2v^_ (22) 

The above example can be generalized to stretched exponentials: For 
f{P) of the form f{P) ~ e~'^^ one obtains after a short calculation 



E 



1/(5-1) 



and 

p(i?)~i?(2-5)/(2<5~2)gaE^/(^-)^ (24) 

where a is some factor depending on 6 and c. In this case the super- 
statistical complex system exhibits stretched exponential tails. 



7 



5 Anomalous diffusion in superstatis- 
tical systems 

We now illustrate that superstatistical systems can exhibit normal as 
well as anomalous transport. This depends on the dynamical proper- 
ties of the model considered. 

For simplicity, we restrict ourselves to a 1-dimensional model. Let 
us again consider locally a 1-dimensional Brownian particle of mass 
m and a Langevin equation of the form 

i, = -jv-^aL{t), (25) 

where v denotes the velocity of the particle, and L(t) is normalized 
Gaussian white noise with the following expectations: 

{L{t)) = (26) 
{L{t)L{t')) = 6{t-t'). (27) 

We assume that the parameters a and 7 are constant for a sufficiently 
long time scale T, and then change to new values, either by an explicit 
time dependence, or by a change of the environment through which the 
Brownian particle moves. Formal identification with local equilibrium 
states in the spatial cells where /3 is approximately constant (ordinary 
statistical mechanics at temperature yields during the time scale 
T the relation [35] 

iv') = f = l^ (28) 
27 pm 

or 

f3=-^. (29) 

Again, we emphasize that after the time scale T, 7 and a will take on 
new values in a stochastic way. During a time interval of the order of 
T, the probability density P{v^ t) obeys the Fokker-Planck equation 

dP divP) 1 od'^P 



with the local stationary solution 



P(.v\P) = \l^e^pi-ll3mv'}. (31) 
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In the adiabatic approximation, valid for large T, one asumes that the 
local equilibrium state is reached very fast so that relaxation processes 
can be neglected. Within a cell in local equilibrium the correlation 
function is given by [35] 

C{t-m = {v{tHt')) = ^e-^\'-'r (32) 

It is now interesting to see that the long-term invariant distribution 
P{v), given by 

/>oo 

P{v) = / f{P)P{v\P)dp (33) 
Jo 

depends only on the probability distribution of /3 = and not on 
that of the single quantities 7 and cr^. This means, one can obtain 
the same stationary distribution (j33p from different dynamical models 
based on a Langevin equation with fluctuating parameters. Either 7 
may fluctuate, and cj^ is constant, or the other way round. On the 
other hand, the superstatistical correlation function 

poo 1 poo 

C{t -t')= f{l3)C{t - t'\f3)d(i = - / fm-^e-^\'-''\df3 (34) 

JQ ^ Jo 

can distinguish between these two cases. The study of correlation 
functions thus yields more information for any superstatistical model. 

Let us illustrate this with a simple example. Assume that a fluc- 
tuates and 7 is constant and that P = -^-^ is ^^-distributed. Since 7 
is constant, we can move the exponential e"''''*"*'' out of the integral 
in eq. (j34p . meaning that the superstatistical correlation function still 
decays in an exponential way: 

C(t-0 ~e-^'*"*''- (35) 

On the other hand, if a is constant and 7 fluctuates and /? is still 
X^-distributed with degree n, we get a completely different answer. In 
this case, in the adiabatic approximation, the integration over P yields 
a power-law decay of C(t — t'): 

C{t-t') r^\t-t'\-i, (36) 

where 

r? = ^-l (37) 
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Note that this decay rate is different from the asymptotic power law 
decay rate of the invariant density P{v), which, using ([3T]) and ([33]) . 
is given by P{v) ~ with 

1 n 1 , ^ 

+ (38) 



q-1 2 2 
Now let us proceed to the position 



rt 

x{t) = I v{t')dt' (39) 







of the test particle. One has 



{x\t)) = f f {v{t')v{t"))dt'dt". (40) 
Jo Jo 

Asymptotic power-law velocity correlations with an exponent r/ < 1 
are expected to imply asymptotically anomalous diffusion of the form 

{X\t))r^e (41) 

with 

a = 2-r]. (42) 

This relation simply results from the two time integrations. 

It is interesting to compare our superstatistical model with other 
dynamical models generating anomalous diffusion. Plastino and Plastino 
and Tsallis and Bukmann[371 study a generalized Fokker-Planck equa- 
tion of the form 

= -^{F{x)P{x,t)) + D^P{x,tr (43) 

with a linear force F{x) = ki — k2X and 7^ 1. Basically this model 
means that the diffusion constant becomes dependent on the proba- 
bility density P. The probability densities generated by eq. (03]) are 
(/-Gaussians with the exponent 

q = 2-u. (44) 

The model generates anomalous diffusion with a = 2/(3 — q). Assum- 
ing the validity of a = 2 — r}, i.e. the generation of anomalous diffusion 
by slowly decaying velocity correlations with exponent 57, one obtains 

V = ^. (45) 
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On the other hand, for the x^'Superstatistical Langevin model one 
obtains by combining eq. (j37p and (|38p the different relation 

.= 2^. (46) 

Interesting enough, there is a distinguished g- value where both models 
yield the same answer: 

q = 1.453 ^fi = 'q = 0.707 (47) 

These values of q and r/ correspond to realistic, experimentally ob- 
served numbers, for example in defect turbulence (see section 9). 



6 From time series to superstatistics 

We now want to be more practically orientated and apply supersta- 
tistical techniques to some complex systems (of whatever kind) where 
we do not know the equations of motion, and neither the distribution 
/(/3), but do have some information in form of a measured time series. 
Suppose an experimentally measured scalar time series u{t) is given. 
Our goal is to test the hypothesis that it is due to a superstatistics and 
if yes, to extract /(/3). First we have to determine the superstatistical 
time scale T. For this we divide the time series into N equal time 
intervals of size Ai. The total length of the signal is tmax = NAt. We 
then define a function K{At) by 

Here (• • ■)toAt — si Ito~^^^ ■ ■ ■ dt denotes an average over an interval 
of length At starting at to- The integration result fluctuates for each 
value of to and is averaged by the integral over to- u denotes the 
average of u. Now assume the simplest case, that our complex system 
dynamics basically arises out of a superposition of local Gaussians on 
some unknown time scale T. How can we extract T? We should be 
looking for the special value At = T where 

k{T) = 3. (49) 

Clearly this condition defining the superstatistical time scale T simply 
reflects the fact that we are looking for locally Gaussian behaviour in 
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Figure 1: Determination of the superstatistical time scale T from the inter- 
section with the hne k = 3 for turbulent Taylor-Couette flow, 6 = 2\j = 
0, 1, 2, . . . , 7 (from top to bottom). 



the time series, which implies a local flatness of 3. If At is so small that 
only one constant value of u is observed in this interval, then of course 
K,{A.t) = 1. On the other hand, if At is so large that it includes the 
entire time series, then we obtain the flatness of the distribution of the 
entire signal, which is larger than 3, since superstatistical distributions 
are generically fat-tailed. Inbetween, there should be a distinguished 
time scale where k = 3. 

Fig. 1 shows the function K(At) for an example of a time se- 
ries that has been studied in [2], the longitudinal velocity difference 
u{t) = v{t + 5) — v{t) \n Bt turbulent Taylor-Couette flow on a scale 5. 
For each scale 5 the relevant superstatistical time scale T leading to lo- 
cally Gaussian behaviour can be extracted as the intersection with the 
line K = 3. These time scales T are to be compared with the relaxation 
times of the dynamics, which can be estimated from the short-time 
exponential decay of the correlation function C{t — t')) = {u{t)u{t')) . 
One obtains the result that the ratio T/j~^ is pretty large and in- 
creases with Reynolds number [2]. This time scale separation is in- 
deed the deeper reason why superstatistical models of turbulence work 
quite weU [EIEIIII]. 

Next, given a general signal u{t) we are interested in the analysis of 
the slowly varying stochastic process P{t). Since the variance of local 
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'lognormal' 
'chi-square' 




0.01 0.1 1 10 

Figure 2: Probability distribution as extracted from the measured tur- 
bulent time series in [2]. 



Gaussians ~ e 2^" is given by /3 ^, we can determine the process 
P{t) from the time series as 

= ^ , ,2 (50) 

We can then easily make a histogram of /3(to) for ah values of to, thus 
obtaining the probability density 

Fig. 2 shows this probability density for our example of turbulent 
time series. Motivated by our consideration in section 3, the data are 
compared with a x'^-distribution, inverse x^-distribution and lognor- 
mal distribution, all having the same mean and variance — 
as the experimental data. Clearly the lognormal distribution yields the 
best fit. Indeed, the cascade picture of energy dissipation in turbu- 
lent flow strongly suggests that lognormal superstatistics should be 
relevant, with P being related to a suitable power of the fluctuating 
energy dissipation rate \27\ However, again let us mention that 
other complex systems can generate a completely different type of 
superstatistics. Still the same methods apply. 

For superstatistics to be a good approximation we need the variable 
pit) to change very slowly as compared to u{t). This is indeed the 
case for our turbulence example, as can be seen in Fig. 3. 
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Figure 3: Typical evolution of u{t) (solid line) and (dashed line) for the 
turbulent Taylor-Couette flow. 



7 Overview of applications 

The superstatistics concept has been apphed to many different com- 
plex systems. Let us first give a short overview of applications, after 
that we will treat four examples in more detail. Rizzo and Rapisarda[161 
[25] study experimental data of wind velocities at Florence airport and 
find that superstatistics does a good job. Swinney et al. [2l|38] extract 
lognormal superstatistics for turbulent flow between counterrotating 
disks. Paczuski et al.p!3] study data of solar flares on various time 
scales and embedd this into a superstatistical model based on x^- 
superstatistics. Human behaviour when sending off print jobs might 
also stand in connection to such a superstatistics [39] . Bodenschatz et 
al.[40l HU |42] and Pinton et al. [43l [44j have detailed experimental 
data on the accelerations of single test particles in a turbulent flow, 
which are well described by lognormal superstatistics [9| 1101 [TT]. The 
statistics of cosmic rays is well described by x^-superstatistics, with 
n = 3 due to the three spatial dimensions |15j. In mathematical finance 
superstatistical techniques are well known and come under the heading 
'volatility fluctuations', see e.g.[T71 [181 119]. Possible applications also 
include granular media, which could be described by different types 
of superstatistics, depending on the boundary conditions j45j- The 
observed fat tails of solar wind speed fluctuations |46] could also be 
related to a superstatistical model. Hydroclimatic fluctuations have 
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been analysed using the superstatistics concept (2^. Briggs et al. [H] 
apply a superstatistical model to observed train delays on the British 
rail network. On the theoretical side, Chavanis[6] points out analo- 
gies between superstatistics and the theory of violent relaxation for 
collisionless stellar systems. Abul-Magd [20] applies superstatistics to 
random matrix theory. Luczka and Zaborek[47) have studied a sim- 
ple model of dichotomous fluctuations of f3 on different time scales 
where everything can be calculated analytically. Mathai and Haubold 
[48j investigate a link between superstatistics and fractional reaction 
equations. 



We now treat some examples of applications in more detail. We start 
with the recent Lagrangian turbulence applications [11] . Over the past 
few years there has been experimental progress [401 ffll H2| [H] in 
tracking single test particles advected by a turbulent flow. This area 
of research is called 'Lagrangian turbulence'. To theoretically model 
Lagrangian turbulence, one may first start from a Gaussian turbu- 
lence model, the Sawford model [49^ I50j . This model considers the 
joint stochastic process {a{t),v{t),x{t)) of an arbitrary component of 
acceleration, velocity and position of a Lagrangian test particle embed- 
ded in the turbulent flow, and assumes that they obey the stochastic 
differential equation 



8 Lagrangian turbulence 
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e: average energy dissipation 

Co,ao: Lagrangian structure function constants 

al variance of the velocity distribution 

R\ = \fvZ Taylor scale Reynolds number. 
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For our purposes it is sufficient to consider tlie limit — > oo, 
wliich is a good approximation for large Reynolds numbers. In that 
limit the Sawford model reduces to just a linear Langevin equation for 
the acceleration 



a = -7a + aL{t) (54) 

with 

7 = (55) 

^3/2 

a = ^u-^/^e. (56) 
2ao 

Note that this is a Langevin equation for the acceleration, so the 
meaning of the variables is slightly different as compared to the case of 
an ordinary Brownian particle, where the Langevin equation describes 
the velocity. In practice, the acceleration is measured as a velocity 
difference on a very small time scale. 

Unfortunately, the Sawford model predicts Gaussian stationary 
distributions for a and v, and is thus at variance with the recent mea- 
surements |4H HH] , which provide evidence for distributions with fat 
tails. So how can we extend the Sawford model to make it physically 
realistic? 

As said before, the idea is to generalize the Sawford model with 
constant parameters to a superstatistical Sawford model with fluctuat- 
ing ones. To construct a superstatistical extension of Sawford model, 
one replaces the constant energy dissipation e by a fluctuating one. 
It is assumed to be lognormally distributed. Moreover, one extends 
the model to include all 3 components of the velocity and accelera- 
tion, as well as contributions from a fluctuating enstrophy (rotational 
energy) surrounding the test particle. From this new theory [TT| ex- 
cellent agreement with the experimental data is obtained, see Fig. 4 
for an example. One obtains not only the correct 1-point probabil- 
ity distributions, but also good agreement for the decay of correlation 
functions, the observed statistical dependencies between acceleration 
components, and scaling exponents. 



9 Defect turbulence 

Let us now consider another physically relevant example, so-called 'de- 
fect turbulence'. Defect turbulence shares with ordinary turbulence 
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Figure 4: Predicted and measured probability density of a component of the 
small-scale velocity difference of a Lagrangian test particle in a turbulent flow. 
The dashed line is based on lognormal superstatistics [llj. The experimental 
data are from |43j . 



only the name, otherwise it is very different. It is a phenomenon 
related to convection and has nothing to do with fully developed hy- 
drodynamic turbulence. Consider a Raleigh-Benard convection ex- 
periment: A liquid is heated from below and cooled from above. For 
large enough temperature differences, interesting convection patterns 
start to evolve. An inclined layer convection experiment [51^ [521 |T3] 
is a kind of Raleigh-Benard experiment where the apparatus is tilted 
by an angle (say 30 degrees), moreover the liquid is confined between 
two very narrow plates. For large temperature differences, the convec- 
tion rolls evolve chaotically. Of particular interest are the defects in 
this pattern, i.e. points where two convection rolls merge into one (see 
Fig. 5). These defects behave very much like particles. They have a 
well-defined position and velocity, they are created and annihilated in 
pairs, and one can even formally attribute a 'charge' to them: There 
are positive and negative defects, as indicated by the black and white 
boxes in Fig. 5. 

The probability density of defect velocities has been quite precisely 
measured jl4J. As shown in Fig. 6, it quite precisely coincides with 
a g-Gaussian with q ~ 1.46. The defects are also observed to exhibit 
anomalous diffusion. Their position X(t) roughly obeys an anomalous 
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Figure 5: Convection rolls and defects (black and white boxes) as observed 
in the experiment of Daniels et al. [14j. 
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Figure 6: Measured probability density of defect velocities and fit with a 
g-Gaussian with g = 1.46. 
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e = 0.08 e = 0.17 




Figure 7: Measured anomalous diffusion of defects for two different values 
of the non-dimensional temperature differences e between upper and lower 
plate. 

diffusion law of tfie type 

{X\t))^t-, (57) 

where 1.33 (see Fig. 7). 

The simple superstatistical model of section 5 with fluctuating ef- 
fective friction 7 makes sense as a very simple model for the defect 
velocity v. While ordinary Brownian particles have constant damping 
due to Stokes' law 7 = where u is the kinematic viscosity of the 

liquid, p is its density, m is the mass of the particle and a is the radius 
of the particle, defects are no ordinary particles: They are nonlinear 
excitations and have neither a well-defined mass m nor a well-defined 
radius a. Thus one expects that there is an ensemble of damping 
constants 7 which depend on the topology of the defect and its fluc- 
tuating environment. In particular, the fastest velocities result from 
circumstances in which the defect is moving in a local environment 
with only a very small effective local damping 7 acting. The driving 
forces L{t) are hardly damped during such a time interval, and lead to 
very large velocities for a limited amount of time, until another region 
with another 7 is reached. The result are g-Gaussians, as shown in 
Fig.6, and anomalous diffusion. In good approximation this system is 
described by x^-superstatistics. 
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10 Statistics of cosmic rays 

Our third example is from high energy physics. We will proceed to 
extremely high temperatures, where (similar as in defect turbulence) 
particles are created and annihilated in pairs. We will apply supersta- 
tistical techniques to high energy collision processes on astrophysical 
scales, leading to the creation of cosmic ray particles that are ul- 
timately observed on the earth. The idea to apply superstatistical 
techniques to the measured cosmic ray spectrum was first presented 
in [15], based on some earlier work in [53]. 

Experimental data of the measured cosmic ray energy spectrum 
are shown in Fig. 8. Also shown is a curve that corresponds to a 
prediction of a superstatistical model. Up to energies of 10^^ eV, the 
measured flux rate of cosmic ray particles with a given energy E is 
well fitted by a distribution of the form 

^'^'-^■ (i+ifa-'iW/..-.. - (^«> 

E is the energy of the particles, 



E = a/ c^p^ + c^p^ + c^pI + m^c^, (59) 
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b = (kf)^^ is an effective inverse temperature variable, and C is a 
constant representing the total flux rate. For relativistic particles the 
rest mass m can be neglected and one has E ~ c\p\. The distribution 
(j58p is a ^-generalized relativistic Maxwell-Boltzmann distribution in 
the formalism of nonextensive statistical mechanics |31) . The factor 
E'^ takes into account the available phase space volume. As seen in 
Fig. 8, the cosmic ray spectrum is very well fitted by the distribution 
()58p if the entropic index q is chosen as q = 1.215 and if the effective 
temperature parameter is given by kT = b^^ = 107 MeV. 

The above effective temperature is of the same order of magnitude 
as the so-called Hagedorn temperature Th [SH |55] , an effective tem- 
perature well known from collider experiments. The Hagedorn tem- 
perature is much smaller than the center-of-mass energy Eq'ms of a 
typical collision process and represents a kind of 'boiling temperature' 
of nuclear matter at the confinement phase transition. It is a kind of 
maximum temperature that can be reached in a collision experiment. 
Even largest Eqms cannot produce a larger average temperature than 
Th due to the fact that the number of possible particle states grows 
exponentially. 

Let us now work out the assumption that the power law of the 
measured cosmic ray spectrum is due to fiuctuations of temperature. 
Assume that locally, in the creation process of some cosmic ray parti- 
cle, some value of the fluctuating inverse temperature P is given. We 
then expect the momentum of a randomly picked particle in this re- 
gion to be distributed according to the relativistic Maxwell-Boltzmann 
distribution 

Pirn = zip-f'''^''- 

Here p{E\P) denotes the conditional probability of E given some value 
of p. We neglect the rest mass m so that E = c\p\ . The normalization 
constant is given by 

poo n 

= E^e-f^'^dE = ^. (61) 

Now assume that /3 is x'^-distributed. The observed cosmic ray dis- 
tribution at the earth does not contain any information on the local 
temperature at which the various particles were produced. Hence we 
have to average over all possible fluctuating temperatures, obtaining 
the measured energy spectrum as the marginal distribution 

/•oo 

p{E) = / p{E\(3)f{f3)d(3. (62) 
Jo 
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The integral (f62]) with /(/3) given by ([6]) and p{E\P) given by ([60|) is 
easily evaluated and one obtains eq. (j58p with 



2 



(63) 



9 = 1 + 



n + 6 



and 



6 



/9o 



(64) 



4-3g 



where /3o is the average inverse temperature. 

The variables Xi in section 3 describe the independent degrees of 
freedom contributing to the fluctuating temperature. At very large 
center of mass energies, due to the uncertainty relation, the probed 
volume is very small, and all relevant degrees of freedom in this 
small volume are basically represented by the 3 spatial dimensions into 
which heat can flow, leading to a fluctuating effective temperature in 
each creation process of cosmic ray particles. The 3 spatial degrees of 
freedom yield n = 3 or, according to eq. ([63ll . 



For cosmic rays Ecms is very large, hence we expect a g-value that is 
close to this asymptotic value. The fit in Fig. 8 in fact uses q = 1.215, 
which agrees with the predicted value in eq. ([65]) to about 3 digits. 



Our final example leaves the area of classical physics and turns to 
a more practical problem that almost everybody has experienced in 
the past. Trains, buses, planes, etc. are often delayed! A statistical 
analysis of train delay data in the UK was recently performed in [21]. 
One observes probability densities of delays that are g-exponentials, 
i.e. that have power-law tails (see Fig. 9 for an example). 

We start with a very simple theoretical model for the local depar- 
ture statistics of trains. The waiting time distribution until departure 
takes place is simply given by that of a Poisson process [35] 



Here t is the time delay from the scheduled departure time, and f3 
is a positive parameter. The symbol P(t|/3) denotes the conditional 
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(65) 



11 Statistics of train delays 



P{t\p) = (3e-P\ 



(66) 
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Figure 9: Delay statistics of train departures at Swindon station. The trains 
are heading for London Paddington. 

probability density to observe the delay t provided the parameter /? 
has a certain given value. Clearly, the above probability density is 
normalized. Large values of (5 mean that most trains depart very well 
in time, whereas small /3 describe a situation where long delays are 
rather frequent. 

The above simple exponential model becomes superstatistical by 
making the parameter [5 a fluctuating random variable as well. These 
fluctuations describe large-scale temporal variations of the British rail 
network environment. For example, during the start of the holiday 
season, when there is many passengers, we expect that [3 is smaller 
than usual for a while, resulting in frequent delays. Similarly, if there 
is a problem with the track or if bad weather conditions exist, we 
also expect smaller values of (3 on average. The value of (3 is also be 
influenced by extreme events such as derailments, industrial action, 
terror alerts, etc. 

The observed long-term distribution of train delays is then a mix- 
ture of exponential distributions where the parameter (3 fluctuates. 
If [3 is distributed with probability density and fluctuates on a 

large time scale, then one obtains the marginal distributions of train 
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delays as 

p{t) = / f{(3)pmd(3 = / f{P)(3e-^'. (67) 
Jo Jo 

Again, a x^-distribution of n degrees of freeedom makes sense for 
P, leading to g'-exponential waiting time distributions of the form 

p{t) ^ {1 + b{q - l)t)^^ (68) 

where q = l + 2/(n + 2) and b = 2/?o/(2 — q). Our model generates 
g-exponential distributions of train delays by a simple mechanism, 
namely a x^-distributed parameter P of the local Poisson process. 
This is in good agreement with the recorded delay data of the British 
rail network p4| . 

Typical (/-values obtained from our fits for various stations are in 
the region (7 = 1.15... 1.35 (see Table 1). Hence 

n = - 2 (69) 

is in the region 4 ... 11. This means the number of degrees of freedom 
influencing the value of /? is just of the order we expected it to be: A 
few large-scale phenomena such as weather, seasonal effects, passenger 
fluctuations, signal failures, repairs of track, etc. seem to be relevant. 
In general, it makes sense to compare stations with the same q (the 
same number of external degrees of freedom of the network environ- 
ment): The larger the value of b, the better the performance of this 
station under the given environmental conditions. 
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station 




Q 




b 


code 


Bath Spa 


1 


.195 


0, 


.209 


BTH 


Birmingham 


1 


.257 


0, 


.271 


BHM 


Cambridge 


1 


.270 


0, 


.396 


CBG 


Canterbury East 


1 


.298 


0, 


.400 


CBE 


Canterbury West 


1 


.267 


0, 


.402 


CBW 


City Thameshnk 


1 


.124 


0, 


.277 


CTK 


Colchester 


1 


.222 


0, 


.272 


COL 


Coventry 


1 


.291 


0, 


.330 


COV 


Doncaster 


1 


.289 


0, 


.332 


DON 


Edinburgh 


1 


.228 


0, 


.401 


EDB 


Ely 


1 


.316 


0, 


.393 


ELY 


Ipswich 


1, 


.291 


0, 


.333 


IPS 


Leeds 


1 


.247 


0, 


.273 


LDS 


Leicester 


1 


.231 


0, 


.337 


LEI 


Manchester Piccadilly 


1 


.231 


0, 


.332 


MAN 


Newcastle 


1 


.378 


0, 


.330 


NCL 


Nottingham 


1 


.166 


0, 


.209 


NOT 


Oxford 


1 


.046 





.141 


OXF 


Peterborough 


1, 


.232 


0, 


.201 


PBO 


Reading 


1 


.251 


0, 


.268 


RDG 


Sheffield 


1 


.316 


0, 


.335 


SHE 


Swindon 


1 


.226 


0, 


.253 


SWI 


York 


1 


.311 


0, 


.259 


YRK 



Table 1: The fitted parameters q and b for the departure statistics of 23 UK 
stations. 
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